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ABSTRACT 


are  created  in  a  supersonic  mole¬ 
cular  jet  and  are  studied  by  two-color  tine-of-flight  mass  spectroscopy.  The 
clusters’  s"-*-  S  1  intermolecular  vibronic  structures  are  characterized  by  calcu- 

i  0 

lational  modeling  of  the  clusters'  intermolecular  motion.  The  calculations 
include:  1)  an  intermolecular  normal  coordinate  analysis  (NCA)  which  treats  all 
six  van  der  Vaals  modes  under  a  harmonic  oscillator  assumption;  and  2)  a  three 
dimensional  hindered  rigid  rotor  analysis  (3D-HRRA)  which  treats  the  intermole¬ 
cular  torsional  motion.  Agreement  between  calculation  and  experiment  is  excel¬ 
lent  for  binding  energy,  symmetries,  and  van  der  Yaala  mode  energies.  The  cluster 
spectra  and  calculated  intermolecular  modes  are  compared  to  those  of  bensene  (CH*) 
reported  previously.  A  major  conclusion  of  this  work  is  that  the  clusters  behave 
rigidly  with  regard  to  Internal  rotation  of  the  cluster  subunits  and  that  the 
clusters  possess  unique  equilibrium  geometries.  The  Internal  torsional  motion 
is  oscillatory  and  that  motion  is  constrained  by  an  orlentationally  dependent 
Intermolecular  potential  whose  barrier  height  is  of  the  order  of  the  cluster 
binding  energy. 
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Clusters  of  benzene  (CD^)^  and  (CF^)^ 
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I.  Introduction* 

van  der  Waals  (vdW)  complexes  of  aromatic  molecules  with  hydrocarbon 
solvents  form  a  class  of  suoramolecular  systems  whose  intermolecular  bonding 
has  several  interesting  features.  First,  the  intermolecular  interaction  is 
small  and  is  dominated  by  long-range  dispersive  attractions  and  short-range 
exchange  repulsions.  These  two  features  allow  the  interaction  to  be  modeled 
using  an  intermolecular  potential  of  known  functional  form.  For  example,  the 
interaction  can  he  modeled  by  additive  atom-atom  potentials  set  in  a  Lennard- 
Jones  or  an  Exponential-Six  format.1  Second,  the  interaction  potential  surface 
only  changes  slightly,  if  at  all,  upon  electronic  excitation  of  the  cluster 
chromophore.  Thus,  small  spectral  shifts  of  the  chromophore  electronic  trans¬ 
ition  and  weak  intensity  intermolecular  vdV  mode  vlbronic  transitions  are 
observed.2'^  Third,  the  interaction  results  in  the  formation  of  specific 
minimum  energy  cluster  configurations.  These  geometries  are  Interesting  since 
they  give  insight  into  the  nucleatlon  processes  and  solvation  geometry  occurring 
in  both  gas  and  condensed  phase  systems.1  Fourth,  the  small  binding  energy  of 
the  complex  and  the  low  frequency  vdV  vibrational  modes  are  important  since  they 
play  an  essential  role  in  intramolecular  vibrational  redistribution  (IVR)  and 
vibrational  predissociation  (VP)  cluster  dynamic  energy  transfer  processes.4  And 
fifth,  cluster  structure,  binding  energy,  vdV  modes  and  dynamics  (IVR  and  VP) 
are  eesentially  dependent  on  the  actual  cluster  structure  and  the  intricacies 
of  the  intermolecular  interaction.^'^ 

In  this  paper,  we  report  the  spectroscopic  results  of  the  benzene  deu- 
terated  methane  (benCCD^)  and  benzene  carbon  tetrafluoride  (benfCF^),) 
clusters  together  with  calculated  modeling  of  selected  cluster  charactenaf:’?. 
The  spectroscopic  results  include  the  w  *■  w  vlbronic  spectra  of  the  clusters 
m  their  respective  benzene  constituent  6^  regions. 
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The  calculated  results  include  the  geometry,  the  binding  energy,  and 
the  full  eigenvalue/eigenvector  intermolecular  vibrational  structure  for  the 
electronic  ground  state  of  each  cluster.  The  intermolecular  ground  state  vi¬ 
brational  structure  is  modeled  by  two  methods:  (1)  an  intermolecular  normal 

coordinate  analysis  (NCA)  which  determines  all  six  intermolecular  vdW  mode 
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fundamentals  under  a  harmonic  oscillator  assumption;'  and  (2)  a  three  dimen¬ 
sional  hindered  rigid  rotor  analysis  (3D-HRRA)  for  which  an  anisotropic 
perturbation  field  is  applied  to  the  cluster  solvents  considered  to  be  rigid 
internal  rotors  in  the  resulting  potential  well.  The  intermolecular  torsional 
mode  structure  is  thereby  obtained  under  a  hindered  rigid  rotor  approximation. 

The  calculated  ground  state  vibrational  structure  for  each  cluater  is 
compared  to  its  experimental  vlbronic  spectrum  and  to  the  structure  observed 
and  calculated  for  the  bensene  methane  (benCCH^)  cluster.1’^  Within  this 
framework,  vlbronic  selection  rules  which  govern  the  clusters'  S1  ♦  SQ  transi¬ 
tions  are  discussed  along  with  the  ramiflcationa  of  Hersberg-Teller  (H-T)  coupling 
for  the  observation  of  "forbidden"  vlbronic  transitions. 

The  motivation  for  studying  these  systems  concerns  the  elucidation  of 
the  vdV  torsional  structure.  Two  limiting  cases  can  be  proposed  in  regard  to 
the  torsional  structure  of  ben(CH^)j,  ben(CD4)1t  and  benlCF^.  In  one 
case,  the  clusters  possess  free  internal  rotation  between  the  cluster  solute 
and  solvent.  The  cluster  solvent  (CH^,  CD^,  or  CF^)  rotates  freely  in  three 
dimensions  against  the  bensene  framework  and  the  system  is  considered  nonrigid. 

In  the  other  case,  the  clusters  possess  torsional  oscillations  for  which  the 
cluster  solvent  librates  against  the  bensene  frame  with  a  residence  time  lone 
enough  to  give  rise  to  "vibration  like"  motion  in  a  rigid  molecule  regime. 

Elucidation  of  the  vdW  torsions  is  of  particular  interest  in  these 
systems  since  the  aforementioned  limiting  cases  pose  questions  regarding  the 
actual  physics  governing  the  torsional  structure:  (1)  do  the  clusters  possess 
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free  internal  rotation  between  the  cluster  solute  and  solvent  or  do  they  possess 
torsional  oscillations  for  which  the  lowest  eigenstates  behave  nore  or  less  har¬ 
monically;  and  (2)  in  either  case,  what  is  the  dependence  of  the  internolecular 
potential  upon  the  relative  orientation  of  the  cluster  constituents? 

II.  Experimental  Procedures. 

The  apparatus  and  techniques  employed  in  this  study  are  largely  the  same 
as  those  previously  reported.1’^  The  ben(CD^)1  and  ben(C?4)1  spectra  are 
recorded  using  a  pulsed  supersonic  molecular  jet  expansion  in  combinstion  with 
two-color  time-of-flight  mass  spectroscopy  (2-color  TOFMS).  Two  Nd*^/YAG  pumped 
dye  lssers  are  used  to  probe  the  cluster  spectra.  An  LBS  698  dye  laser  whose 
output  is  frequency  doubled  and  mixed  with  the  Hd+^/YAG  1 .064  micron  fundamental 
is  used  to  probe  the  clusters'  S.  *■  sA  transition  in  the  isolated  bensene  61 
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region  ('b^  ♦  'a^;  *  •*•*)•  Subsequent  ionisation  of  the  clusters  is  accom¬ 
plished  using  a  R590  dye  laser  whose  output  is  frequency  doubled  and  mixed  with 
the  lfd*^/YAG  1.064  micron  fundamental.  The  ionisation  laser  output  is  set  at 
45112  cm-1.  A  mixture  of  deuterated  methane  or  carbon  tetrafluorlde  in  helium 
is  placed  in  line  with  liquid  bensene  in  a  trap  at  room  temperature.  The  three 
component  mixture  is  expanded  using  a  pulsed  molecular  jet  nossle  having  a  500 
micron  orifice  while  maintaining  100  pslg  backing  pressure.  Apparatus  chamber 
pressure  is  maintained  at  or  below  5  x  10”®  torr  during  the  experiments. 

III.  Theoretical  Considerations. 

The  RCAs  are  conducted  using  the  same  methods  as  described  in  our  pre¬ 
vious  publication  on  vdV  cluster  vibronic  structure.^  For  ben(CD^)1t  the 
calculated  cluster  ground  state  geometry,  force  field,  and  binding  energy  are 
taken  as  those  of  ben(CH. Only  the  masses  are  changed  for  the  deuterium/ 
hydrogen  substitutions  in  the  NCA.  For  benCCF^)^  the  cluster  ground  sta*e 
geometry  and  binding  energy  are  calculated  via  intermolecular  energy  minimisa¬ 
tion  employing  the  methods  previously  described.1'^  The  NCA  is  conducted  usir.r 


the  modeled  geometry  and  its  corresponding  force  fielu.  The  potential  para¬ 
meters  for  the  C-F  Ucp  =  6.08  x  10?  cm'1 /A12;  Ccp  =  1.1 3  x  105  cm'1 /A6) 

and  H-F  (A  *  7.77  x  106  cm_1/A;  C,™  =  3.12  x  104  cm'VA)  atom-atom  interactions 
HF 

are  determined  from  parameters  and  methods  described  elsewhere. 

The  3D-HRRA  centers  upon  analyzing  the  effects  of  the  intermolecular 
potential  on  the  cluster  solvent  rotational  degrees  of  freedom.  In  the  absence 
of  a  perturbing  field  from  the  cluster  solute,  the  cluster  solvent  rotates  freely 
in  three  dimensions  against  a  fixed  cluster  solute  framework.  Application  of  a 
specified  perturbation  field  results  in  the  eventual  restriction  of  this  free 
rotation  to  torsional  oscillation  as  the  perturbation  field  magnitude  ia  increased 
that  is,  rotation  ceases  and  vibrational  oscillatory  motion  begins  as  the  resi¬ 
dence  times  of  the  solvent  in  the  torsional  potential  well  becomes  longer. 

In  the  3D-HRRA,  both  cluster  constituents  are  assumed  internally  rigid, 
their  respective  internal  geometries  remaining  constant  and  at  "equilibrium”. 

The  cluster  solvents  are  taken  as  the  rotating  portions  of  the  clusters  since 
their  zero  field  rotational  constants  are  orders  of  magnitude  greater  than  those 
of  the  cluster  solute.  The  hindered  solvent  rotations  are  presupposed  to  con¬ 
tribute  to  the  spectral  features  observed  along  with  the  intermolecular  bending 
and  stretching  modes. 

The  3D-HBRA  involves  setting  up  a  molecule  fixed  coordinate  system  (x, 
y,z)  and  a  space  fixed  coordinate  system  (  (  •  n  , c  )  on  the  cluster  solvent  as 
shown  in  Pigure  1.  Both  systems  have  their  origins  at  the  nuclear  center-of- 
mass  of  the  cluster  solvent.  The  molecule  fixed  coordinate  system  is  chosen 
such  that  its  principle  axes  lie  along  the  rotational  axes  of  the  solvent 
tetrahedron  (T^  point  group  symmetry).  The  space  fixed  coordinate  system  is 
chosen  to  be  at  the  Euler  angles  9  •  43  degrees,  $  *  0  degrees,  and  x *  180 
degrees  with  respect  to  the  molecule  fixed  coordinate  system.  The  relative 


orientation  of  these  two  systems  ia  chosen  so  that  the  orientational  dependence 
of  the  perturbation  field  can  be  expressed  as  a  simple  function  of  the  curvi¬ 
linear  coordinates  9  ,  $  ,  and  x • 

An  analytical  expression  for  the  potential  governing  the  torsional  notion 
can  be  derived  by  analyzing  the  solvent  rotations  in  each  of  the  rotational  coor¬ 
dinates.  Since  a  solvent  rotation  is  periodic  over  360  degrees,  the  potential  can 
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be  expressed  by  an  even  function  Fourier  series  written  as 

V  » 

V(q)  m  +  2  A  cos(p q)  (1) 

P-1  P 

with  q  a  function  of  the  rotational  coordinates  and  V  taken  as  the  barrier 

o 

height  to  internal  rotation.  The  function  q  is  determined  by  the  potential  dimen¬ 
sionality  and  the  number  of  minima  in  the  potential  along  each  of  the  rotational 
coordinates.  The  first  cosine  term  in  the  series  describes  the  major  form  of  the 
potential,  and  the  rest  of  the  terms  in  the  series  "fine  tune"  the  potential 

shape.  Since  the  coefficients  A  in  these  terms  are  assumed  to  be  small,  the 

y 

series  can  be  truncated  at  its  first  cosine  term  without  sacrificing  the  general 
potential  shape.  Doing  this  yields  the  mathematically  convenient  potential  form: 

V 

V(q)  -  y  d  *  cos (<?)).  (2) 

The  symmetry  of  the  cluster  solvents  can  be  employed  to  determine  the 
dimensionality  of  and  the  number  of  minima  in  the  potential  and,  hence,  the 
function  q .  Rotating  the  cluster  solvents  against  the  benzene  frame  by  90 
degrees  in  0  or  x  ,  270  degrees  in  x  »  or  180  degrees  in  the  $>  coordinate 
results  in  equivalent  maxima.  Equivalent  minima  occur  at  0  degrees  rotation 
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in  9,  ♦,  and  x  and  at  180  degrees  rotation  in  9  and  x*  The  maxima  have 
the  same  energy  relative  to  the  equivalent  minima  in  all  three  rotational 
coordinates.  Generalizing  these  relative  displacements  to  all  coordinate 
values,  one  finds  that  the  potential  magnitude  depends  upon  the  relative 
displacements  of  all  three  rotational  coordinates  as, 

q  -  29  +  2X  +  (3) 

Inserting  eq.  (3)  into  eq.  (2)  yields  the  three  dimensional  potential  function 
used  in  the  calculations: 

V 

V(  9,  <j>,  X)  *  "j[l  -  cos(29  ♦  2x  +  ♦)]  (4) 

in  which  Vq  is  the  harrier  height  to  internal  rotation.  The  rotational  wave- 
functions  chosen  as  the  basis  set  for  the  calculation  are  the  rigid  rotor  sym¬ 
metric  top  wavefunctions  which  depend  on  the  curvilinear  coordinates  and  on  the 

quantum  numbers  J,  k  and  m.  Under  zero  field  conditions,  these  wavefunctions 
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are  solutions  to  the  spherical  top  Schrodinger  equation, 

ljkm>  -  E^otljkm>  (5) 

with  eigenvalues  corresponding  to  those  of  the  spherical  top  cluster  solvent. 
Upon  application  of  the  perturbation  field,  the  eigenvalues  are  obtained  via 
diagonalization  of  the  energy  matrix  possessing  elements  of  the  form 

<j’k  m’|  'ryJ2  ♦  V(  9  ,  ♦  ,  x  )|  J  *  m  ), 

■n  d 


(6) 


using  eq.  (l)  for  V(  Q,  $  ,  x)  with  Vq  as  an  adjustable  parameter.  In  the 
calculations  performed,  the  only  possible  nonzero  off-diagonal  matrix  elements 
occur  with  the  following  restrictions: 


J  -  J  | 12 

(7) 

k  - k  |=2 

(8) 

•  ll 

m  -  a  =  1 . 

(9) 

Diagonal  matrix  elements  contain  only  the  zero  field  spherical  top  energies 
and  a  V  J2  potential  term. 

The  resulting  matrix  is  diagonalized  for  a  selected  value  of  Vq 
yielding  eigenvalues  corresponding  to  the  solvent  torsional  eigenstates  at  the 
specified  perturbation.  The  perturbation  is  varied  until  a  reasonable  fit  with 
the  experimental  spectrum  is  obtained.  The  calculations  are  performed  on  a 
Cyber  205  computer  using  a  basis  set  consisting  of  680  wavefunctions  to  ensure 
convergence  of  the  lowest  eigenstates  at  their  proper  eigenvalues.  Only  the 
torsional  structures  of  ben(CH^)^  and  ben(CD^)1  are  calculated  since  their 
respective  rotational  constants  are  large  enough  that  the  coupling  of  the  rota¬ 
tional  levels  does  not  require  an  extremely  large  basis  set  to  insure  converg¬ 
ence.  Matrix  elements  are  determined  via  numerical  integration  using  a  non- 
adaptive  integration  routine.  The  matrix  is  prediagonalized  into  a  tridiagonal 

Q 

form  using  orthogonal  similarity  transformations  and  diagonalized  using  an 
implicit  QL  method.1*^  The  rotational  constants  used  for  methane  and  deuterated 
methane  are  5*2  cm"1 ,  and  2.6  cm-1 ,  respectively. 

TV.  Results. 

A.  Ben(CD ).. 

4  i 

Fig.  2  and  Table  1  present  the  ben(CD^)1  vdW  cluster  spectrum  in 
the  region  between  38552.4  cm"1  and  38678.6  cm"1  using  2-color  TOFMS.  The 
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cluster  61  transition  is  "red  shifted"  by  41 .2  cm  *  with  resoect  to  the 
o 

benzene  61 .  The  bathychromic  shift  i3  comparable  to  that  of  ben(CH.)1  (41  cm-1) 

indicating  that  the  binding  energy  difference  between  S1  and  So  for  ben(CD4)1 

is  nearly  identical  to  that  of  benCCH^.  Ten  pronounced  internolecular 

vibronic  features  are  observed  to  the  "blue"  of  the  cluster  61  .  As  in  the 

o 

case  of  the  ben(CH^)1  cluster,  no  features  are  observed  in  the  symmetry  for¬ 
bidden  benzene  0°  region.  Thus,  the  cluster  geometry  must  possess  at  least  a 
three-fold  axis  of  symmetry. 

The  calculated  ben(CD^)^  geometry  (fig.  3)  is  assumed  to  be  the  same 
as  that  calculated  for  ben(CH^)^.  The  geometry  possesses  C3v  point  group 
symmetry.  In  this  geometry,  the  CD^  center-of-mase  lies  3*47  A  above  the 
benzene  molecular  plane  along  the  three- fold  rotational  axis.  The  ground  state 
binding  energy  in  this  configuration  is  calculated  at  540  cm”1 .  Using  the 
cluster  "red  shift"  of  41.2  cm”1,  the  excited  state  binding  energy  is  calcu¬ 
lated  to  be  581  cm”1 . 

The  NCA  reveals  six  vdW  vibrations  (fig.  3  and  Table  1),  two  being 

two-fold  degenerate.  The  ground  state  vibrational  energies  are  75  cm  for 

the  vdW  stretch  s  (a,),  15  cm-1  for  the  bends  b  (e),  and  20  cm-1  (t  (a~)) 
z  i  xy  z  z 

and  64  cm  1  for  the  vdW  torsions.  The  eigenvector  normal  modes  (fig.  3) 

transform  in  an  identical  fashion  to  those  calculated  for  bentCH^.  The  vdW 
stretching  mode  transforms  as  the  translation  of  the  cluster  constituents  away 
frcm  one  another  along  the  z  (3-fold)  axis.  The  vdW  bends  transform  as  some 
combination  of  cluster  constituent  translations  in  opposite  directions  perpendi¬ 
cular  to  the  3-fold  axis  in  the  xy  plane.  One  vdW  torsion  transforms  as  a  rot¬ 
ation  of  the  cluster  constituents  about  the  3-fold  axis  in  opposite  directions. 
The  remaining  two  vdW  torsions  transform  as  rotations  about  orthogonal  axes 
perpendicular  to  the  3-fold  axis. 

The  3D-HRRA  results  are  also  included  in  Table  1  for  the  lowest  eigen¬ 
states  in  the  torsional  manifold.  The  torsional  mode  structure  is  calculated 


for  V  =  "500  cm-1  and  B  s  2.6  cm-1.  This  perturbation  results  in  a  reasonable 
o 

fit  for  the  torsional  features  observed  experimentally.  The  torsional  "zero 
point  energy"  is  75  cm"1 .  Two  distinct  torsional  manifolds  result  from  the 
calculations:  one  manifold  has  eigenvalues  grouped  quartically  (nearly  four 
fold  degenerate  located  at  approximately  21  cm  1 ,  46  cm  1 ,  and  75  cm  1  above 
the  zero  point  energy)  and  the  other  manifold  has  eigenvalues  grouped  octally 
(nearly  eight  fold  degenerate  located  at  approximately  68  cm-1  above  the 
torsional  zero  point  energy). 

For  comparison,  the  observed  vibronic  features  for  ben(CH^)^  are  repro¬ 
duced  in  Table  1  along  with  the  results  of  the  NCA  and  the  3D-HRRA.  The  NCA 
results  for  ben(CH^)^  are  those  reported  previously.  The  3D-HRRA  torsional 

structure  is  calculated  for  V  *  300  cm"1  and  B  M  5*2  cm"1.  This  perturbation 

o 

is  chosen  since  it  is  assumed  that  the  barrier  to  internal  rotation  is  nearly 
identical  in  the  two  systems  as  they  only  differ  by  isotopic  substitution  and 
have  the  same  electronic  structure.  Additionally,  using  the  same  potential 
barrier  of  electronic  origin  for  both  CH^  and  CD^  clusters  provides  a  check  on 
the  validity  of  the  model  in  predicting  the  torsional  mode  structure  of  the 
clusters.  For  ben(CH^)^  the  torsional  "zero  point  energy"  is  96  cm  1 .  As  in  the 
case  of  the  ben(CD^)^  calculations,  two  distinct  torsional  manifolds  result:  one 
manifold  with  eigenvalues  grouped  quartically  at  about  29  cm"1 ,  54  cm-1  and  73  cm 
above  the  torsional  zero  point  energy  and  the  other  manifold  with  eigenvalues 
grouped  octally  at  about  89  cm"1  above  the  zero  point. 

B.  Ben(C7j)j. 

Fig.  4  and  Table  2  present  the  ben(CF^)^  cluster  spectrum  recorded 
using  2-color  TOFMS  in  the  region  between  38578.6  cm"1  and  38702.6  cm”1.  Un¬ 
like  the  ben(CDA)1  and  ben(CH^)1  cluster  6^  transitions,  the  ben(CF^)1  cluster 
6^  is  blue  shifted  by  6.1  cm"1  with  respect  to  the  benzene  6\  The  small  hypso- 
chromic  shift  indicates  that  the  binding  energies  in  the  Sq  and  states  are 


nearly  identical  with  the  ground  state  binding  energy  being  slightly  greater. 

The  relative  displacement  between  the  two  potential  surfaces  is  also  small  since 


only  5  vdW  transitions  are  observed  and  their  intensities  decrease  abruptly  at 

about  50  cm-1  above  the  cluster  6^  origin.  Mo  cluster  spectrum  is  observed  in 

the  forbidden  benzene  0°  region  indicating  that  the  cluster  possesses  at  least 

o 

a  three-fold  rotation  axis.  This  result  is  not  unexpected  since  the  same  obser¬ 
vation  is  made  for  the  ben(CH^)^  and  ben(CD^)^  systems. 

The  ground  state  configuration  and  vdW  eigenvalues/eigenvectors  are 
shown  in  fig.  5  and  Table  2.  Only  one  minimum  energy  geometry  is  calculated 
for  the  cluster.  The  geometry  has  point  group  symmetry  with  the  CF^ 
center-of-mass  at  3.43  A  above  the  benzene  molecular  plane  along  the  3-fold 
axis.  The  cluster  ground  state  binding  energy  is  1064  cm-1. 

The  NCA  reveals  six  vdW  vibrations.  Their  ground  state  vibrational 
energies  are  69  cm'1  for  the  vdW  stretch  s  (a.),  11  cm-1  for  the  vdW  bends 
b  (e),  and  13  cm'1  (*  (a2))  and  36  cm'1  (t  (e))  for  the  vdW  torsions. 

The  eigenvector  normal  modes  transform  as  the  translational  and  rotational  repre¬ 
sentations  of  the  C^v  point  group  as  indicated  and  in  the  same  manner  as  those 
for  ben(CH4)1  and  ben(CD^)j. 

V.  Discussion. 

In  comparing  the  calculated  ground  state  vdW  vibrational  structure  and 

experimental  vdW  vibronic  structure,  we  assume  that  the  intermolecular  potential 

surfaces  of  the  clusters  studied  are  identical  in  both  the  5.  and  S  electronic 

1  o 

states.  Additionally,  we  utilize  the  group  theoretical  arguments  developed  in 
our  previous  publication^  on  vdW  cluster  vibronic  structure  to  assign  and  under¬ 
stand  the  observed  cluster  spectra. 


A.  Ben(CD^)^. 


The  calculated  ground  state  vdVf  vibrations  (NCA)  of  benCCD^  and 
the  experimental  vibronic  spectrum  is  shown  in  fig.  2  and  Table  1.  As  previously, 
the  assignments  are  made  by  direct  comparison  between  the  calculations  and  the 
experimental  vibronic  spectra.  The  intense  feature  at  22.5  cm”1  to  the  blue  of 
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the  cluster  6^  corresponds  to  the  t2  torsion  calculated  at  20  cm”1.  Thus  this 
feature  and  its  observed  overtones  at  48.4-  cm  1  ,  70.5  cm  1 ,  and  100  cm  1  are 

assigned  to  6^  vdW  combination  bands  6gtz^a2^0’  ^otz^a2^o’  ^otz^a2^o’  and 

6 1 t  (a~)4*  The  features  at  32  cm"1  and  62.7  cm'1  to  the  blue  of  the  cluster 
o  z  2  o 

6^  correspond  to  the  first  and  third  overtones  of  the  2-fold  degenerate  vdW 

bending  modes  whose  fundamental  is  calculated  at  15  cm  1 .  They  are  assigned  to 

the  61  combination  bands  61b_^>(e)^t  and  6^b  (e)4,  respectively.  The  feature 

at  58.1  cm-1  to  the  blue  of  the  cluster  6^  is  assigned  to  a  2-fold  degenerate 

t  vdW  torsions/cluster  61  combination  band,  61 t  (e)1.  The  t  „  torsion  fund- 
xy  o  o  xy  o  xy 

amental  is  calculated  at  64  cm”1 .  With  this  assignment  and  that  of  the  t 

z 

torsion,  the  feature  at  84  cm”1  is  assigned  to  the  61 1  (e^t  (a.,)1  combina- 

o  xy  o  z  do 

tion  band.  The  vdW  stretch  is  calculated  at  75  cm”1.  This  mode  is  identified 

in  the  cluster  vibronic  spectrum  at  72.9  cm”1  to  the  blue  of  the  cluster  6^ 

origin  and  is  thus  assigned  as  Under  this  scheme,  we  assign  the 

features  at  93*3  cm”1  and  97*5  cm”1  as  61s  (a.^b _ (e)1  and  61s„(a1  J1 1 

o  z  1  o  xy  o  ozioz 

(a-)1  combination  bands. 

2  O 

The  torsional  structure  resulting  from  the  3U-HRRA  for  ben(CD^)^,  (Table 
1),  confirms  the  torsional  assignments  made  using  the  NCA.  The  quartically  grouped 
torsional  levels  at  21  cm”1 ,  46  cm”1  and  75  cm”1  correspond  to  those  associated 
with  the  tz  torsion  and  its  overtones  in  the  cluster  spectrum  at  22.5  cm  1 ,  48.4 
cm”1  and  70.5  cm”1.  The  octally  grouped  torsional  levels  at  68  cm”1  correspond 
to  the  t  torsion  assigned  at  58.1  cm  1 . 

The  correspondence  between  these  levels  and  those  calculated  using  the 
NCA  can  be  understood  as  follows.  In  the  3D-HRRA,  four  symmetrically  equivalent 
minima  exist  in  the  torsional  potential  surface.  These  minima  correspond  to  the 
four  ways  of  placing  the  solvent  tetrahedron  upon  the  solute  with  a  tetrahedral 
face  toward  the  solute  molecular  plane.  If  the  barrier  between  these  minima  i3 
infinitely  high,  penetration  of  the  local  wavefunctions  through  the  barrier 


separating  the  potential  minima  does  not  take  place*  In  this  case,  each  well 
contains  eigenstates  corresponding  to  intermolecular  torsions  which  occur  with 
small  amplitude  about  each  potential  well  minimum.  Since  the  four  potential 
wells"  are  identical  in  shape  and  depth,  a  four-fold  "structural"  degeneracy 
exists  in  which  all  four  potential  wells  contain  identical  torsional  structure. 
Thus,  for  example,  a  nondegenerate  torsional  eigenstate  actually  has  a  four¬ 
fold  structural  degeneracy,  etc. 

Considering  the  actual  situation  in  which  the  potential  barrier  is  fi¬ 
nite,  tunneling  occurs  and  the  "structural"  degeneracy  is  lifted  via  interaction 
of  the  local  wavefunctions  through  the  potential  barrier.  The  eigenstate  split¬ 
ting  due  to  this  tunneling  may,  or  may  not,  be  observed  depending  on  the  experi¬ 
mental  resolution  and  the  relative  difference  between  the  eigenstate  energy  and 
the  barrier  height.  Splitting  of  the  "structural"  degeneracy  in  the  lower  por¬ 
tion  of  the  potential  well  is  minimal  unless  the  barrier  is  low.  Based  on  the 
calculations,  energy  level  splittings  of  the  first  few  sets  of  eigenstates  should 
not  be  observed  unless  the  barrier  is  below  150  cm”1. 

At  moderate  barriers  (about  300  cm”1),  the  eigenstates  in  the  lower 
portion  of  the  well  are  nearly  degenerate  and  behave  more  or  less  harmonically. 
For  all  practical  purposes,  we  can  assume  these  levels  to  be  degenerate.  The 
torsional  level  structure  can  then  be  determined  from  the  eigenstates  in  one  of 
the  minima.  Thus,  the  3D-HRRA  torsional  structure  calculation  simplifies  into 
the  NCA.  Physically,  the  cluster  can  be  considered  at  least  "semi-rigid"  in  the 
respect  that  it  has  a  definable  equilibrium  configuration.  The  potential  energy 
barrier  separating  one  minimum  from  the  others  in  the  potential  surface  is  lar<a 
and  may  be  of  the  order  of  the  cluster  binding  energy. 

The  vibronic  structure  in  both  the  ben(CH4)1  and  ben(CD4)1  spectra 
suggests  that  the  systems  are  more  or  less  rigid.  None  of  the  observed  vibrcr.i? 
state  energies  follow  a  free  rotor  formalism  for  which  the  energy  level  structure 


ig  described  by  eq.  (2).  If  the  systems  behaved  nonngidly,  the  free  rotor  eigen¬ 
states  would  lie  at  approximately  2.6 J(J*l)  cm  and  5 .2J(J*i>  cm  for  ber.  "T,  , 
and  ben(CH^) , ,  respectively,  and  a  lJm  ^  selection  rule  would  govern  the  Tar.si- 
tions. 

In  the  ben(CH.),  and  ben(CD,).  systems,  the  torsional  mode  shifts 

4  1  4  1  Z 

by  17.6  %  upon  deuterium  substitution  of  the  cluster  solvent.  Both  the  NCA  and 

3D-HRRA  models  predict  a  28.6%  frequency  shift  for  t  ,  while  the  free  rotor 

z 

model  would  predict  a  50%  frequency  shift  for  t  upon  deuteration.  Considering 
that  only  one  mode  is  taken  into  account  and  the  mode  couplings  may  be  different 
in  the  protonated  and  deuterated  clusters,  we  conclude  that  torsional  tunneling 
does  not  take  place  in  the  lower  levels  of  the  bentCH^)^  and  ben(CD^)j  cluster 
potential  wells  and  that  the  clusters  can  be  considered  to  be  rigid. 

The  negligible  isotopic  red  shift  in  the  b^  bends  also  suggests  that 
the  cluster  is  rigid  and  that  the  entire  observed  spectrum  is  not  solely  due  to 
internal  rotation.  The  same  rationale  holds  for  the  vdW  stretch:  its  observa¬ 
tion  also  dispels  the  notion  that  only  rotor  modes  occur  in  the  spectrum.  No 
experimental  isotopic  shift  can  be  determined  for  the  stretching  mode,  however, 
since  it  is  not  observed  in  the  ben(CH^)1  system.  Theoretically,  the  mode  should 
red  shift  by  8.5%  upon  deuteration  of  the  cluster  solvent.  The  ben(CH4)1  stretch 
should  then  he  at  about  80  cm-1  in  the  vibronic  spectrum  based  upon  the  observe  i 
stretch  in  the  benCCD^  system.  (The  ben(CH^)1  stretch  is  calculated  by  the 
NCA  to  be  at  82  cm-1)* 


The  ben(CD^)^  cluster  spectrum  is  also  substantially  richer  than  the 

ben(CH^)1  spectrum  in  the  respect  that  both  the  vdW  stretch  and  t  torsions 

are  observed  along  with  well  developed  t  torsion  and  bend  progressions.  T1-: 

z 

is  probably  due  to  more  favorable  Franck-Condon  factors  in  the  benCCD^)^  case 


Av.'V'f 
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*9  -.9  'he  case  for  ben(CH.),,  H-T  coupling  influences  the  vibromc 

4  T 

intensities  (selection  rules)  in  the  ben(CD4)1  spectrum.  The  nontotally 

symmetric  modes  t  ,  t  and  k  all  appear  in  the  spectrum  wi'h  lv  -  _*i  ,  *_2. 
z  xy  xy 

♦t,  ...  selection  rules  as  can  be  seen  in  fig.  2  and  Table  1. 

B.  BeeiCF^) | • 

Comparison  of  the  calculated  ground  state  vdW  vibrations  (NCA)  of 

ben(C?4)j  and  the  experimental  cluater  spectrum  is  presented  in  fig.  4  and 

Table  II.  The  intense  feature  at  17.4  cm-1  to  the  blue  of  the  ben(CF. )  61 

4  1  o 

corresponds  to  the  tg  torsion  calculated  at  13  cm-1.  Thus,  this  feature  and  its 

first  ovartone  at  35*4  cm-1  are  aaaigned  to  the  61  vdV  combination  banda  61 

o  o 

t  (a-)1  and  6^t  (a,,)2.  The  shoulder  feature  at  15.7  cm”1  to  the  blue  of 

S  c  O  O  S  c  0 

the  cluster  6^  corresponds  to  the  vdV  2-fold  degenerate  b^  bend  fundamental 

calculated  at  11  cm'1.  Taking  this  feature  aa  the  combination  band  6^  (e)1 

o  xy  o 

1  2  *1 

leads  to  the  assignment  of  the  b  bend  first  overtone  6  b__ (e)  at  32.3  cm  . 

xy  o  xy  o 

Finally,  the  spectral  feature  at  39.5  cm”1  in  the  ben(CF^)t  spectrum  corre¬ 
sponds  to  the  t  torsion  calculated  at  36  cm”1 •  Thus,  the  feature  la  assigned 
xy 

to  the  combination  band  6*t  (e)1. 

o  xy  o 

The  ben(CF4)^  spectrum,  like  the  ben(CD4)1  and  ben(CH4)1  spectra, 
suggests  that  the  system  is  rigid.  Since  the  rotational  constant  for  CF^  is 
small  (ca.  0.18  cm”1),  the  free  rotor  energy  level  structure  should  appear  at 
about  .37  cm”1  intervals.  This  structure  is  not  observed.  Instead,  the  spec¬ 
trum  possesses  oscillatory  torsional  structure  commensurate  with  the  NCA  theore* 
ical  predictions.  These  theoretical  and  experimental  results  demonstrate  the 
rigidity  of  the  systems. 

The  S.  *■  S  excitation  of  the  ben(CF.),  cluster  involves  very  little 
change  in  cluster  geometry.  In  this  spectrum,  the  progression  intensities  ie- 
crease  dramatically  at  energies  greater  than  the  6^  cluster  origin.  In  the 
lowest  two  vdW  mode  progressions,  the  intensities  decrease  approximately 
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between  -he  -'de  f  indanentals  and  the  first  overtones.  Further-.ore,  the  nigh 
„r.«r, ry  vd i  stretch  calculated  at  *3  cm"'  is  not  observed.  At  "hese  high  e-er- 
gies,  the  Franck-Condon  factors  must  be  very  small. 

■"he  observation  of  the  nontotally  symmetric  vdV  b  bends  and  the  t^ 

*  Jr  — 

torsion  progressions  with  a  ^v  *  ' ,  *2  ...  selection  rule  suggests  that  mter- 

electronlc  state  mixing  ( H-T  coupling)  is  an  important  contributor  to  the  mode 

intensity  mechanism.  Aa  in  the  ben(CH^)1  and  ben(CD4)1  cases,  H-T  coupling 

becomes  apparent  in  the  low  lying  vdV  modes  of  the  bentCF^  system.  In  fact, 

the  interelectronic  state  mixing  is  substantial  in  the  ben(CF4)1  system  and  can 

be  demonstrated  by  the  observation  of  the  well  defined  band  fundamental  at  15.7 

cm"1  and  by  the  observation  of  the  t  torsion  fundamental  (39.5  cm"1)*  Both 

r7 

of  these  latter  features  should  not  be  obsarvsd  if  H-T  coupling  is  not  present. 
YI.  3s— ary  and  Comclnmloma. 

Two-color  TOFMS  and  supersonic  molecular  jet  techniques  have  been 
employed  to  study  the  Sj  <-  SQ  vibronic  spectra  of  ben(CD4)t  and  benfCF^),  vdk 
clusters.  These  studies  reveal  detailed  Information  regarding  the  geometries, 
the  intermolecular  energetics,  and  the  physical  nature  of  the  vdW  interactions 
present  in  the  systems.  The  experimental  observations  demonstrate  that  the 
clusters  are  at  least  aeml-rigld  systems  possessing  unique  equilibria  geometries 
and  that  the  intermolecular  motion  present  in  the  systems  is  oscillatory.  ~hr^ 
comparison  of  these  spectra  with  those  of  ben(rH^)1  previously  studied  m  this 
laboratory,  the  "isotopic"  shifts  observed  in  the  intermolecular  vibronic  struc¬ 
tures  demonstrate  that  all  the  low-lying  intermolecular  eigenstates  are  nearly 
harmonic.  They  are  not  admixtures  of  oscillatory  intermolecular  bends  and 
stretches  and  free  intermolecular  rotation  as  would  occur  if  the  clusters  be¬ 
haved  r.onrigidly. 
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•nent  of  the  intermolecular  motion.  Both  the  NCA  and  3D-HRRA  models  indicate  that 
the  intermolecular  torsional  motion  is  oscillatory  and  that  the  motion  is  con¬ 
strained  by  an  orientationally  dependent  intermolecular  potential.  Detailed 
analysis  of  the  intermolecular  interaction  reveals  that  the  torsional  motion  is 
governed  by  a  three  dimensional  potential  possessing  torsional  barriers  on  the 
order  of  the  cluster  binding  energy.  At  this  barrier  magnitude,  torsional  tun¬ 
neling  in  the  lover  portions  of  the  potential  is  minimal  and  the  residence  time 
of  the  cluster  constituents  in  a  particular  veil  minimum  is  substantially  longer 
than  the  timescale  of  experimental  observation. 

The  assignment  of  geometry  and  identification  of  the  intermolecular  modes 
using  the  theoretical  models  has  also  proven  useful  in  understanding  the  physics 
governing  the  spectroscopic  properties  of  the  clusters.  The  observation  of  non- 
totally  symmetric  intermolecular  mode  fundamentals  and  combination  bands  in  the 
vlbronic  spectra  suggest  that  interelectronic  state  mixing  (H-T  coupling)  is  an 
important  factor  in  the  overall  intensity  mechanism  governing  the  transitions. 
Additionally,  the  Interelectronic  state  mixing  is  most  important  for  the  lov- 
lylng  Intermolecular  bending  and  torsional  modes. 

Both  the  NCA  and  the  3D-HRRA  adequately  model  the  Intermolecular  modes 
In  the  systems  studied;  hovever,  the  NCA  is  the  more  useful  here  due  to  its 
mathematical  simplicity  and  passive  treatment  of  all  six  intermolecular  modes 
in  a  single  calculation.  Unlike  the  3D-HRRA,  the  NCA  does  not  require  detailed 
—  Priori  knovledge  of  the  potential  surface  along  each  coordinate  of  motion. 
Furthermore,  the  large  amount  of  energy  level  coupling  occurring  in  the  3D-HRF.A 
at  moderate  barriers  suggests  that  the  model  is  best  used  in  systems  possessmi? 
free  or  nearly  free  internal  rotation. 


16 


The  rigidity  of  these  clusters  can  he  qualitatively  understood  by  con¬ 
sidering  the  nature  of  the  solvent  rotational  dynamics.  Any  rotational  state 
for  which  J  f  0  has  a  displacement  in  the  e  coordinate:  5  and  <  displacements 
may  or  may  not  be  present  depending  on  the  values  of  the  k  and  m  quantum  numbers. 
Given  the  geometry  of  the  system  then,  low  barrier  one-dimensional  rotation 
(solvent  rotating  around  the  cluster  three-fold  C-X  bond  axis)  is  not  possible. 
Basically  the  two  dimensional  hindering  in  benCCX^)^  clusters  is  insufficient 
to  produce  the  one  dimensional  rotation  occuring  in  "attached  tops"  such  as 
toluene. 1 1 


We  wish  to  thank  Professor  W.  Klemperer  for  many  helpful  discussions 
and  suggestions  concerning  the  rigidity  of  vdV  molecules. 
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TABLE  I 


vdw  spectral  features  tn  b«n(CDi>1  and  ban(CH4>1  6q  regions 


and  calculated 

ground  scire  vdW 

modes  (refer  to 

fig.  2). 
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c  )  Va  lues  ir 

i  parent neses  represent 

"structural"  degeneracy  calculated  in  1D-HRRA. 

i)  Fro*  references  1  and  ). 


*  Assignment 


TABLE  II 


vdW  spectral  features  in  ben(CF^) ^  6^  region  and  calculated 
ground  state  vdW  modes  (refer  to  fig*  4). 


ben(CF4) 

*  1  a) 

Energy 


Observed _  NCA 


0(38614.7) 
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11 

(V 

17.4 

13 
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32.3 
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a)  Energies  are  reported  in  c*  *  relative  to  the  6*  cluster  origin. 

b)  vdW  node  representations  as  per  fig.  4. 


FICURB  CAPTIORS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Stereoscopic  projection  of  solvent  tetrahedron  showing 
relative  orientation  of  the  molecule-fixed  (x,y,z)  and 
space-fixed  (€,n,C)  coordinate  systems  used  in  the  3D- 
HRRA.  The  (x,y,z)  coordinate  system  is  located  at  0  * 

45°.  $  •  0°,  and  x  *  180°  from  the  (C »n »C )  coordinate 
system. 

Two-color  TOFMS  «-  S  spectrum  and  calculated  ground 
state  vdV  modes  (RCA)  of  ben(CD^)^.  Energy  scale  is  relative 
to  ben(CD.).  61  transition  (38567.4  cm-1).  Feature  posi- 

4  10 

tions  and  assignments  as  per  Table  I  and  fig.  3* 

Calculated  ground  state  minimum  energy  configuration  (a)  and 
RCA  eigenvalue/eigenvector  vdW  modes  (b)  -  (g)  for  ben(CD^)1. 
Cluster  symmetry  is  C^y  with  an  equilibrium  intermolecular 
distance  of  3-47  X.  Eigenvectors  are  normalized  and  displayed 
at  2x  magnification  (2  £  total  displacement). 

Two-color  TOFMS  S^  SQ  spectrum  and  calculated  ground 

state  vdW  modes  (RCA)  of  ben(CF^)^.  Energy  scale  is  relative 

to  ben(CF.)«  61  transition  (38614.7  cm-1).  Feature  positions 
4  1  0 


and  assignments  as  per  Table  II  and  fig.  5. 


Figure  5 


Calculated  ground  state  minimum  energy  configuration  (a)  and 
NCA  eigenvalue/eigenvector  normal  modes  (b)  -  (g)  for  benCCF^. 
Cluster  symmetry  is  C^  with  an  equilibrium  intermolecular 
distance  of  3.43  %.  Eigenvectors  are  normalized  and  displayed 
at  2x  magnification  (2  X  total  displacement). 
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